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Abstract. A large fraction of brown dwarfs and low-mass stars may form by gravitational 
fragmentation of relatively massive (a few 0.1 Mq) and extended (a few hundred AU) discs 
around Sun-like stars. We present an ensemble of radiative hydrodynamic simulations that 
examine the conditions for disc fragmentation. We demonstrate that this model can explain the 
£\J \ low-mass IMF, the brown dwarf desert, and the binary properties of low-mass stars and brown 
^ . dwarfs. Observing discs that are undergoing fragmentation is possible but very improbable, as 
the process of disc fragmentation is short lived (discs fragment within a few thousand years). 
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1. Introduction 

The formation of brown dwarfs (BDs) and low- mass stars is not well understood (Whit- 
worth et al. 2007). Low- mass objects are difficult to form by gravitational fragmentation 
O ,. of unstable gas, as for masses in the BD regime (< 80 Mj, where Mj is the mass of 
Jupiter), a high density (> 10~ 16 g cm~ 3 ) is required for the gas to be Jeans unstable. 
Padoan k, Nordland (2004) and Hennebelle & Chabrier (2008) suggest that these high 
density cores can be formed by colliding flows in a turbulent magnetic medium. How- 
ever, this model requires a large amount of turbulence, and has difficulty in explaining 
the binary properties of BDs. Additionally, the large number of brown-dwarf mass cores 
that the theory predicts have not been observed. Another way to reach the high densi- 
ties required for the formation of BDs and low-mass stars is in gravitationally unstable 
discs (Whitworth & Stamatellos 2006; Stamatellos et al. 2007b; Stamatellos & Whit- 
worth 2008, 2009a, b). These discs form around newly born stars and grow quickly in 
. mass by accreting material from the infalling envelope. They become unstable if the 
mass accreted onto them cannot efficiently redistribute its angular momentum outwards 
in order to accrete onto the central star (Attwood et al. 2009). BDs are also thought 
to form as ejected embryos from star forming regions, i.e. as a by-product of the star 
t— ( . formation process (Reipurth & Clarke 2001). In this model BDs form the same way as 
low-mass stars, i.e. in collapsing molecular cores, but shortly after their formation they 
• i-h , are ejected from their parental core and stop accreting any further material. Hence, they 
do not realise their potential to become hydrogen-burning stars. In this paper we focus 
on the mechanism of BD formation by fragmentation of gravitationally unstable discs. 
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2. The fragmentation of gravitationally unstable discs 

Discs can fragment if (i) they are massive enough so that gravity overcomes thermal 
and local centrifugal support (Toomre 1964), and (ii) they can cool fast enough so that the 
energy provided by the collapse of a proto-fragment is radiated away and the growth of 
the proto-fragment continues (Gammie 2001; Rice et al. 2005). Analytical and numerical 
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Figure 1. Left: The mass spectrum of the objects produced by disc fragmentation. Most of 
these are BDs (70%); the rest are low-mass stars. The vertical dotted lines correspond to the 
D-burning limit (~ 13 Mj) and the H-burning limit (~ 80 Mj). The red dashed line refers to 
a low-mass IMF with AN/ AM oc M -0,6 . Right: The BD desert. The semi-major axes of the 
objects (at t — 200 kyr) formed in the disc and remained bound to the central star plotted 
against their mass. The bars indicate the minimum and maximum extent of the orbit. There is 
a lack of BD companions close to central star, but there is a population of BDs loosely bound 
to the central star. The red dots correspond to low-mass binary systems formed in the disc. 

studies have shown that the cooling time must be on the order of the dynamical time 
which happens to be similar to the orbital time. 

Massive (a few times 0.1 Mq), extended (> 100 AU) discs exist (e.g. Enoch et al. 
2009; Jogersen et al. 2005; 2009; Greaves et al. 2008; Andrews et al. 2009; Duchene 
this volume). We have performed an ensemble of radiative hydrodynamical simulations 
of such discs. We assume a star-disc system in which the central star has initial mass 
Mj = 0.7Mq. Initially the disc has mass M n = 0.7Mq, inner radius R m = 40 AU, outer 
radius R OVT — 400 AU, surface density £(-R) oc i?~ 7 / 4 , temperature T(R) oc i? -1 / 2 , and 
hence approximately uniform initial Toomre parameter Q ~ 0.9. Thus, the disc is at the 
outset marginally gravitationally unstable. The only parameter that is different between 
different runs is the random noise from which the gravitational instabilities grow in the 
disc. The evolution of the disc is initially followed using SPH, until ~ 70% of the disc mass 
has been accreted, either onto the stars condensing out of the disc, or onto the central 
star; this typically happens within 20 kyr. Then the residual gas is ignored and the long 
term dynamical evolution of the system is followed up to 200 kyr, using an N-body code. 
The energy equation and associated radiative transfer are treated with the method of 
Stamatellos et al. (2007a). The intrinsic radiation of the central star is taken into account 
but its accretion luminosity is ignored (Stamatellos & Whitworth 2009a). The radius of 
the sink representing the star is set to 1 AU. The gravitational instabilities develop quickly 
and the disc fragments within a few thousand years. Typically 5-10 objects form in each 
disc, most of them BDs. The final status of these objects is determined by subsequent 
accretion of material onto them and by their mutual interactions. 

3. The statistical properties of objects formed by disc fragmentation 

The mass spectrum. Most of the objects (~ 70%) formed in the discs are BDs, 
including a few planetary-mass BDs. The rest are low-mass hydrogen-burning stars. The 
typical mass of an object produced is ~ 20 — 30 Mj. The shape of the initial mass 
function (Fig. [TJ left) is similar to what observations suggest, i.e. it is consistent with 
AN/ AM oc M~ a , where a w 0.6 (e.g. in Pleiades a = 0.6 ± 0.11 - Moraux et al. 2003 - 
and in a Orionis a = 0.6 ± 0.1 - Lodieu et al. 2009). 
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The brown dwarf desert. Most of the fragments form in the discs with initial masses 
as low as 3 Mj (Stamatellos & Whitworth 2009c). After a fragment forms it accretes 
mass from the disc, and interacts with the disc through drag forces, and dynamically 
with other fragments. As a result some fragments migrate close to the central star. This 
region is rich in gas, hence these accreting proto-fragments eventually become low-mass 
hydrogen-burning stars. Fragments that form farther out also accrete material from the 
disc but not quite as much; these become BDs. If any of the BDs happen to migrate in 
the region close to the central star they tend to be ejected back into the outer region 
through 3-body interactions. Hence, the disc fragmentation model produces a lack of BD 
close companions to Sun- like stars, i.e the BD desert (Fig.[l] right). The BD desert may 
extend out to 300 AU but it is less dry outside 100 AU. We predict a population of BDs 
at distances from 20 to 5000 AU from the central star. 

Free-floating planetary-mass objects. The disc fragmentation model provides an 
explanation of the existence of free-floating planetary-mass objects (Lucas & Roche 2000; 
Zapatero Osorio et al. 2000). In our model these objects form in the disc by gravitational 
fragmentation and quickly after their formation they are ejected from the system due to 
3-body interactions. Hence, they stop accreting and their mass remains low. 

Low-mass binary statistics. Close and wide BD-BD binaries are common outcomes 
of disc fragmentation. The components of the binary form independently in the disc and 
then pair up. The simulations produce all kinds of low-mass binaries: star-star, star-BD, 
BD-BD, and BD-planetary mass object binaries. Our model predicts a low-mass binary 
fraction of 16%. This is comparable with the low-mass binary fraction observed in star- 
forming regions (e.g. in Taurus : > 20%, Kraus et al. 2006, in the field 15 ± 5%, Gizis 
et al. 2003). We also predict that close low- mass binaries should outnumber wide ones, 
and this seems to be what is observed (Burgasser et al. 2007). Most of the low-mass 
binaries (55%) have components with similar masses (q > 0.7), in agreement with the 
observed properties of low-mass binaries (e.g. Burgasser et al. 2007). The model also 
predicts that BDs that are companions to Sun-like stars are more likely to be in binaries 
(binary frequency 25%) than BDs in the field (frequency 5 — 8%). This trend is consistent 
with what is observed (Burgasser et al. 2005; Faherty et al. 2009). 



4. Numerical Issues 

We use 1.5 x 10 5 particles to represent each disc, which means that the minimum 
resolvable mass (corresponding to the number of neighbours used, i.e. 50 SPH particles) 
is ?» 0.25 M-,. We have also performed simulations using 2.5 x 10 5 and 4 x 10 5 particles. 
In these simulations the growth of gravitational instabilities, and the properties of the 
proto-fragments formed as a result of these instabilities, follow the same patterns as in 
the simulation with lower resolution, but the details of the final outcomes (i.e. the number 
of objects formed and their exact formation positions) are different. 

The minimum Jeans mass (M, = (47r 5 / 2 /24)(cl/(G 3 p) 1 ^ 2 )) in a typical simulation 
using 1.5 x 10 5 particles is M JEANS MIN m 2 Mj. Thus, according to the Bate & Burkert 
(1997) condition, this mass is adequately resolved as it corresponds to ~ 8 x N NEIGH ( a 
minimum factor of 2 x AT NEIGH is recommended). The Toomre mass (Mt — 7rCj/G 2 S) is 
also adequately resolved. The minimum Toomre mass in our simulation is M TOOMHE MIN fts 
2.5 Mj). This mass corresponds to ~ 10 x -ZV NEIGH ( a factor of 6 x iV NEIGH is recommended; 
Nelson 2006). Finally, the vertical disc structure is resolved by at least ~ 5 smoothing 
lengths, satisfying the Nelson (2006) condition. This is important as the disc mainly cools 
in this direction. 
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5. Observing fragmenting discs 

Typically a few BDs form in each fragmenting disc around a Sun-like star. Hence, it 
may be that only ~ 20% of Sun-like stars need to have large and massive unstable discs 
to produce a large fraction of the observed BDs. Assuming that the lifetime of the Class 
phase is 10 5 yr and that in the disc fragmentation scenario the disc fragments and 
therefore dissipates within 10 4 yr, then the probability of observing a fragmenting disc 
around a Class object is only 10%. Then, considering that only 20% of Sun-like stars 
may have such unstable discs, the probability of observing such discs is only 2%. Hence, 
fragmenting discs should be very difficult to discover (Maury et al. 2009). 

6. Conclusions 

Discs can fragment at distances £ 100 AU (Stamatellos & Whiworth 2009a) from 
the central star to form predominately brown dwarfs, but also low-mass hydrogen burn- 
ing stars and planetary-mass objects. Despite the fact that the model does not include 
magnetic fields, radiative feedback due to the accretion luminosity from newly formed 
protostars, and mechanical feedback (i.e. jets), it can reasonably reproduce (i) the shape 
of the low-mass IMF, (ii) the brown dwarf desert, (iii) the binary properties of low-mass 
objects and (iv) the formation of free-floating planetary mass objects. 
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